Direct High-Power Laser Acceleration of Ions for Medical Applications 
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Theoretical investigations show that linearly and radially polarized multiterawatt and petawatt 
laser beams, focused to subwavelength waist radii, can directly accelerate protons and carbon nuclei, 
over micron-size distances, to the energies required for hadron cancer therapy. Ions accelerated by 
radially polarized lasers have generally a more favorable energy spread than those accelerated by 
linearly polarized lasers of the same intensity. 
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Protons and heavier ions are now being used to treat 
cancer at a number of places around the world Ion 
lithography schemes seem to be heading for practical ap- 
plication [2j and fusion research continues to attract con- 
siderable attention and to gain in importance [3j. In ad- 
dition, considerable effort is being devoted to research 
into the fundamental forces of nature, the initiation of 
nuclear reactions and into schemes to treat radioactive 
waste (ij. In these applications, conventional accelera- 
tors (synchrotrons, cyclotrons and linacs) are employed 
which are large and expensive to build and operate. 

To produce and accelerate ions, current plasma-based 
research focuses on the use of thin foils irradiated by fem- 
tosecond laser pulses of intensity > 10 18 W/cm 2 (Bj. Typ- 
ically a laser pulse is incident on the thin foil giving rise 
to an overdense plasma from which the electrons get ac- 
celerated, form a dense sheath on the opposite side and 
generate a quasistatic electric field of strength in excess 
of 10 12 V/m. This superstrong field accelerates the ions 
to tens of MeV over a distance in the /im range [|| . Re- 
cent work Q has shown that proton beams produced by 
this method of target normal sheath acceleration (TNSA) 
may be improved in energy and beam quality by the use 
of foils less than 1 fim in thickness In earlier exper- 
iments, employing thicker foils, a small fraction of the 
energy got converted to proton energy. Furthermore, the 
protons had energy spreads reaching 100%. 

In hadron therapy [^, Hoi . for example, the ions 
are required to have kinetic energies K = 20 — 250 
MeV/nucleon (H+ and He 2+ ) and K = 85 - 430 
MeV/nucleon (C 6+ and 8+ ). The ions also ought to 
have an energy spread AKf/Kf < 1% so that they may 
be focused on the tumor while sparing the neighboring 
healthy tissue. A beam of rectangular cross section is 
also desirable 



In this Letter we study direct laser acceleration config- 
urations of protons and bare carbon nuclei. The aim is to 
make predictions regarding the optimum conditions that 
would lead to the ion energies of interest to hadron ther- 
apy. A source for the ions may be a dedicated electron 
beam ion trap/source (EBIT/EBIS) [T3| from which they 
can be extracted in a well-defined fully ionized charge 



state or an ensemble of fully stripped ions produced by 
laser-solid interaction. We consider a situation in which 
the source has been tailored on a nanoscale 
study the dynamics of an ensemble of N particles having 
normally distributed kinetic energies, using the single- 
particle relativistic Lorentz-Newton equations. Our cal- 
culations show that laser pulses generated by 0.1 — 10 PW 
laser systems accelerate ions directly to energies in the 
ranges required for hadron therapy, provided the laser 
beams are focused to subwavelength waist radii. The ac- 
celerated ions turn out to have high beam quality and 
rectangular and circular cross sections reflecting symme- 
try of the field and shape of the initial ionic distribution. 

Key to generating the ultrastrong accelerating fields 
is focusing the laser beam to a subwavelength waist ra- 
dius. According to recent experiments, a linearly polar- 
ized beam of wavelength A may be focused to a spot of 
size (0.26A) 2 , where the spot size is the area enclosed by 
a contour at which the beam intensity falls to one half its 
peak value. A radially polarized beam may be focused to 
the substantially smaller spot size of (0.16A) 2 

Upon tight focusing, a linearly polarized laser beam 
develops five field components. The nonvanishing field 
components of a beam polarized along the x axis, propa- 
gating along the z axis, of wavelength A and frequency u>, 
are given in [l5| . using the familiar Gaussian-beam pa- 
rameters (waist radius wq, Rayleigh length z r = ttw^/X 
and diffraction angle e = Wo/z r ). The laser power ex- 
pression may be given, to order e 10 , by 
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where (1q is the permeability of free space, Eqi is the elec- 
tric field amplitude, c is the speed of light in vacuum and 
the subscript I stands for linearly polarized. Note that 
Eqi oc \fP~i and that the leading term in Eqi is inversely 
proportional to wq. 

Dynamics of a particle of mass M and charge Q in 
electric and magnetic fields is governed by the equations 
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FIG. 1: (Color online) Schematic geometry of the acceleration 
scenario. The initial coordinates of the ions are randomly 
distributed in a cylinder which models the interaction regime 
of the ion beam and the laser field. The ejected ions then form 
a beam of circular cross section when the ions are accelerated 
b 
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FIG. 2: (Color online) Coordinates (zf,yj) of 5000 particles 
at the end of their trajectories after interaction with linearly 
polarized light, (a)-(c) are for protons H + and (d)-(f) are 
for carbon nuclei 12 C 6+ . The parameters are: A = 1054 nm, 
wo = 0.5A. The equations have been integrated over a time 
interval equivalent to 1000 laser field cycles, approximately 
amounting to 3.516 ps. For further data see Table [I] 



of motion 
dp 
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in which the energy and momentum of the particle are 
given by £ = jMc 2 and p = jMc/3, respectively, with (3 
its velocity scaled by c, and 7 = (1 — /3 2 ) -1 / 2 its Lorentz 
factor. Integrating Eqs. © numerically one obtains (3 
and, hence, 7/ at a later time tf taken equal to many laser 
field cycles. The final kinetic energy is Kf = (jf — l)Mc 2 . 

To model the interaction region of the ion beam and 
the laser we consider an ensemble of N = 5000 parti- 
cles initially randomly distributed in a cylinder of radius 
R c = 10 nm and length L c = 100 nm oriented along the 
z axis and centered on the origin (see Fig. [I]). In order 
to simulate a realistic ion beam extracted from some ion 
source, the particles will be assumed to possess Gaussian 
distributed random initial kinetic energies with a mean 
value K — 10 keV and a spread AK — 10 eV. Initial 
direction of motion of all particles will be in the xz plane 
and at 9i = 10° relative to the pulse propagation direc- 
tion. 



All E components, together with the v x B force, work 
to accelerate a particle and deflect it to varying degrees 
from its initial direction of motion [l| . A particle makes 
the longest excursion along the x axis due to E x being 
the strongest accelerating component, smaller excursions 
occur due to E y and E z . An ensemble of such particles 
forms a beam of rectangular cross section in the yz plane, 
reflecting symmetry of the laser field. This is demon- 
strated in Fig. [2] for interaction with 0.1, 1 and 10 PW 
laser beams. The ion beam cross section increases with 
increasing laser power (Eqi oc y/Pi). The transverse pro- 
ton beam divergence is larger than for carbon. Energy 
gain by protons is also greater than the gain per nucleon 
of carbon, due to the proton's larger charge to mass ratio. 
Further data not estimable from Fig. [5] are collected in 
Table Q] Note that the energies fall within the domain 
required for hadron therapy and their spread is close to 
what would be suitable for such a purpose. Our calcula- 
tions also show that the relative energy spread tends to 
increase approximately linearly with the volume of the 
initial ionic distribution. However, the input (and out- 
put) particle beam may be collimated using externally 
applied electromagnetic fields (see, e.g. Ref. [16J ) . 

For realization, a short pulse consisting of only a small 
number of field cycles would be needed. This conclusion 
is elucidated by showing, in Fig. [3] for a typical member 
of the ensemble, the particle kinetic energy as a function 
of the number of interaction cycles. The number of ac- 
tual interaction field cycles is small and decreases with 
increasing power, since ions accelerated to higher veloci- 
ties leave the focal region faster. Laser-to-particle energy 
conversion occurs in the form of a few violent impulses. 
Note also that, due to its lower velocity, a carbon ion 
interacts with more field cycles than a proton. 

Estimates may be obtained for the transverse and lon- 
gitudinal emittances ex and ej,, respectively, of the par- 
ticle beam. Using ensemble averages of the exit particle 
coordinates and kinetic energies, one calculates ex < 107r 
mm mrad and £j, < 4x 10~ 5 eV s for protons (case of Fig. 



TABLE I: Data related to Fig. [2] Laser power, average 
particle final x coordinate Xf, average final kinetic energy Ki 
and relative kinetic energy spread AK1/K1 at the end of the 
trajectories of protons, (a)- (c), and carbon nuclei, (d)-(f). The 
ions are interacting with linearly (I) polarized light. 





Power 


Xf 


Ki 


AK t /K, 




[PW] 


W 


[MeV/nucleon] 


[%] 


(a) 


0.1 


89.6 ±0.6 


3.77 ±0.05 


1.3 


(b) 


1 


283.0 ± 1.7 


37.39 ± 0.45 


1.2 


(c) 


10 


750.6 ±41.1 


416.5 ± 24.7 


5.9 


(d) 


0.1 


44.9 ±0.3 


0.94 ±0.01 


1.2 


(e) 


1 


141.3 ±0.9 


9.32 ± 0.11 


1.2 


(f) 


10 


434.3 ± 2.5 


89.9 ±0.9 


1.0 




t/T 



FIG. 3: (Color online) Kinetic energy of a typical particle (out 
of the ensembles of Fig. [2} during interaction with linearly 
polarized light. In the figure, T is one laser field period. 



2(c)). A comparison of these figures with results recently 
reported for ion beams produced from ultrathin solid 
foils irradiated by multiterawatt and petawatt linearly 
polarized laser beams 13|, 13] reveals that a good quality 
proton beam may be obtained by direct laser accelera- 
tion. Estimates for carbon yield < 0.37T mm mrad and 
cl < 9 x 10~ 6 eV s (case of Fig. 2(f)). The remaining 
cases in Fig. 2 have better emittances. 

A Gaussian beam of linear polarization is relatively 
easy to generate in the laboratory, while the same thing 
may not be said about a beam of the radially polarized 
variety. The electric field of a radially polarized focused 
(axicon) laser beam has two components, radial E r and 
axial E z , with propagation along the z axis. However, 
only one azimuthally polarized magnetic held component, 
Bg, exists. E z works efficiently to accelerate the particles, 
while E r and Bg help to limit their diffraction [3, E3| ■ 
Full expressions for the axicon field components are given 
in 20]. The power of the axicon beam, to order e 10 , reads 



Pr = 



' ] E 0r 
' CjJ, 



225 re 

+3o y + -(-2 



(3) 



Note that the amplitude E$ r oc \[P~^. On the other hand, 
when the definition of e is used, one finds that the leading- 
term in Eor is independent of wq. While Eqi has a peak 
value beyond which it falls asymptotically to zero, with 
increasing wo, Eo r increases to an asymptotic value. 

Figs. [HandOare similar to Figs. [5]and[31 respectively, 
albeit for particle acceleration by radially polarized laser 
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FIG. 4: (Color online) Same as in Fig. [2] for particles acceler- 
ated by radially polarized light. For further data see Table ITT1 



fields. From Fig. [4] one sees that the particle beam cross 
section in the xy plane is almost circular, reflecting sym- 
metry of the initial distribution. Recall that the radial 
electric field component E r vanishes on the beam axis and 
its magnitude at a typical particle trajectory is smaller 
than the transverse field in the linearly polarized case. In 
addition, E r plays a confining role during the half-cycle 
in which it is radially inward and tends to cause disper- 
sion during the following half cycle. The fact that most 
of the beam power is concentrated in E z , due to tight fo- 
cusing, is also responsible for the extra energy gain from 
the radially polarized beam compared to that obtained 
from a linearly polarized one. Comparison of Figs. [5] and 
[3] also shows that a particle, since it is confined to dom- 
inantly axial motion, interacts with many more radially 
polarized cycles than it does with linearly polarized ones. 

To assess the beam qualities, estimates show that, for 
protons, tT < 37r mm mrad and cl < 4x 10" 5 eV s (case 
of Fig. 4(c)), while for carbon e<r < 7r mm mrad and 
e_L < 4 x 10" 6 eV s (case of Fig. 4(f)). The remaining 
cases in Fig. 4 have better emittances. These results are 
several orders of magnitude lower than their counterparts 
in conventional accelerators. They are also comparable 
to, and sometimes even better than, the correspondin 
figures for ion beams produced from ultrathin foils [l 



TABLE II: Data related to Fig. H Same as Table |H for ions 
interacting with radially (r) polarized light. 





Power 


z f 




AK r /K r 




[PW] 


[A] 


[MeV/nucleon] 


[%] 


(a) 


0.1 


104.7 ± 0.5 


4.26 ± 0.03 


0.7 


(b) 


1 


430.8 ± 2.3 


45.7 ±0.4 


0.8 


(c) 


10 


3347.5 ± 89.6 


532.8 ± 13.3 


2.5 


(d) 


0.1 


48.5 ±0.2 


1.00 ±0.01 


1 


(e) 


1 


173.7 ±0.8 


10.4 ±0.1 


0.8 


(f) 


10 


876.1 ±5.3 


122.0 ± 1.0 


0.8 



4 









— V 

— 12 c 6+ 






(a) P = 0.1 PW 











40 



60 



NO 



LOO 















(b) P = 1 PW 














FIG. 5: (Color online) Kinetic energy of a typical particle (out 
of the ensembles of Fig. [4]| during interaction with radially 
polarized light. 



17]. 



The average kinetic energies gained from interaction 
with radially polarized laser beams are somewhat bet- 
ter than from interaction with the linearly polarized ones 
(compare Tables U and [TTJ) . These energies, with their 
low spread, cover the domain of application in hadron 
therapy. Note that the initial ensemble volume corre- 
sponds to an ion density rij ~ 10 20 cm -3 . Estimations of 
the energies transferred via space-charge effects suggest 
that those are negligible for the parameters applied. One 
could, in principle, start out with ensemble parameters 



that would lower this density down to rn 



10 



17 



cm 



and still keep the accelerated beam quality within the 
limits of utility in ion therapy. Unfortunately, these den- 
sities are many orders of magnitude higher than what is 
available today from conventional ion sources. Further 
design improvements on such machines need to be made 
before they may be used as ion sources for direct accel- 
eration. An alternative source could possibly be a solid 
target of micro- or nanoscale thickness, perhaps blown off 
by a laser pulse, as in the TNSA mechanism [7J, [8| and 
the laser-piston regime [2l|, followed at the appropriate 
time delay by an accelerating pulse of the type discussed 
in this Letter. 

In conclusion, direct laser acceleration of ions is put 
forward as an appealing alternative for utilization in can- 
cer therapy. Radially polarized beams of the required 
power and intensity are not currently available. Em- 
ploying an existing linearly polarized beam, on the other 
hand, may call for a focusing mechanism to bring the ion 
beam spreading down and render it useful in applications. 



Z.H. acknowledges insightful conversations with Jose 
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